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In this work, we experimentally and theoretically explore voltage controlled oscillations occurring in micro- 
beams of vanadium dioxide. These oscillations are a result of the reversible insulator to metal phase transition 
in vanadium dioxide. Examining the structure of the observed oscillations in detail, we propose a modified 
percolative-avalanche model which allows for voltage-triggering. This model captures the periodicity and 
waveshape of the oscillations as well as several other key features. Importantly, our modeling shows that 
while temperature plays a critical role in the vanadium dioxide phase transition, electrically induced heating 
cannot act as the primary instigator of the oscillations in this configuration. This realization leads us to 
identify electric field as the most likely candidate for driving the phase transition. 



I. INTRODUCTION 

Vanadium Dioxide (VO2) has been a material of pro- 
longed scientific interest, due to the plethora of un- 
usual properties associated with the Insulator to Metal 
phase Transition (IMT) occurring just above room 
temperatur^. The large conductivity change ratio, com- 
bined with an accessible transition temperature and rich 
correlated-electron physics^ ^ has made this an attrac- 
tive compound for many researchers. Much attention has 
historically revolved around controversy over the driving 
physics of the phase transition; particularl y whe ther it 
is a Mott transitiorPlEI or Pierels transitioiPl^. How- 
ever, also of interest is the ab ility of the IMT to happen 
on ultrafast (lOOfs) timescal eJ^^ , and the wide range of 
stimuli which can trigger itt^E^. Along these lines, re- 
cent interest has also shifted from purely academic to in- 
dustrial as well; following proposed applications ranging 
from optical devices^^ and hybrid- metamaterial J^^^^' to 
electronic components'^ and data storage'^"!^. With 
this rise of potential applications comes opportunities for 
new avenues of research and development, but also new 
challenges to satisfy the durability and flexibility that 
real-world devices demancP^^ Understanding the role of 
temperature and structural transitions in various VO2 
phenomenon is key to pushing towards potential of ap- 
plications. 

In thi s man uscript we take an interest in the recently 
reportecP^^ phenomenon of self-sustaining oscillations 
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in VO2 bridges. The widespread prevalence of volt- 
age controlled oscillators in electronics makes this phe- 
nomenon an enticing candidate for devices. It is fairly 
well accepted that these oscillations represent a trigger- 
ing of the Insulator-to-Metal Transition (IMT), followed 
by a reseting Metal-to-Insulator Transition (MIT). How- 
ever, despite headway on controlling such oscillations in 
V02^^ 24^ there is still debate over whether the under- 
lying driving mechanism is thermal or electrostatic. In 
literature, VO2 is most often thermally triggered, and 
yet these oscillations appear to respond to foremost to 
applied voltage across the device. The unavoidable pres- 
ence of joule-heating currents through the 2-terminal de- 
vice during operation, coupled with the observed sensi- 
tivity of the oscillations to device temper at ur^^f^, make 
for a contentious situation. 

In our investigation we first experimentally reproduce 
the oscillations discovered by the authors of ReP^. The 
use of a high-performance oscilloscope in our experi- 
ment gives us access to very fine time-resolution data 
which is useful in our modeling. The details and data 
of our experiment are reported in Section |ll[ Follow- 
ing this, we develop a model which replicates and ex- 
plains the observed waveshape in terms of electrically- 
triggered domains. Our model, reported in Section |III[ 
describes a network of electrically and/or thermally trig- 
gered grains. This model is inspired by several previ- 
ously proposed models^^ While these previous mod- 
els also predict avalanche- like transitions under the right 
conditions^, our our model expands on this framework 
to track time-dependent effects and allow possibility of 
a voltage-triggered phase-transition. Alongside voltage 
triggering, we investigate the role of temperature in the 
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FIG. 1. Schematic of the circuit diagram used to reproduce 
oscillations, depicting the VO2 device as a variable resistance 
and capacitance. Inset shows optical photograph of a sample 
device. The AI2O3 substrate is 330/im thick, and is mounted 
on a glass cover slip. All experiments are performed at room 
temperature. 
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FIG. 2. Applied voltage and voltage across the device as a 
function of time. 



oscillations, and importantly, we find while a voltage- 
driven picture replicates experimental data - thermal 
heating alone is quantitatively and qualitatively unable 
to explain the structure of the observed oscillations. Nev- 
ertheless, device temperature does affect oscillations, and 
thermal co-factors to the voltage triggering are needed 



to reproduce aspects of the data. In Section IV we dis- 
cuss how the percolative transition of VO2 affects the 
shape of the MIT, and what this means regarding ef- 
fective medium within the phase-coexistence region. In 
Section [V| we conclude the manuscript with an overview 
of our results, and an outlook on possible directions for 
VO2 research and application. 



II. VOLTAGE CONTROLLED OSCILLATIONS. 

Our investigation begins with replicating VO2 oscilla- 
tions using the procedure reported by Kim et.aP^^. A 
device consisting of a 10/im x 10/im VO2 bridge between 
two large (^400/im) metal (Ni:Au) electrodes (Figure [T 
inset ) is hooked in series with a limiting resistor {Rext 
and voltage source (Vapp)- This setup is shown schemat- 
ically in Figure [l] Although we do not intentionally add 
external capacitance, the presence of such Cext instru- 
ments and cables is unavoidable - and should be included 
in the effective circuit. The applied voltage Vapp is a 
transient square pulse (between 1/is to 1ms) from a Agi- 
lent function generator (model 33120A) riding on top of 
a constant bias voltage {Vbias = 12^). This is shown as 
the black trace in Figure [2j giving a peak applied voltage 
of 22V. The voltage across the device (V^, as shown by 
the blue trace in Figure |2| is monitored with a LeCroy 
(model wavepro 7-zi) oscilloscope, which allows for high 
time-resolution (40GS/s) sampling resolution even over 
millisecond-long pulses. 

In this configuration, the VO2 device functions es- 
sentially as a capacitor with a variable internal shunt 



resistance Kd- The capacitance Cd is primarily fixed 
by device geometry, although variations of the dielec- 
tric constant of VO2 throughout the phase transition 
(such as have been shown in the context of memory- 
capacitance^^ and VO2 hybrid- metamaterials^^) may 
have small effects - and we discuss this later in Sec- 
tion |IVB| The pre-pulse steady-state starting voltage is 
y D=ybias^D / i^D + ^ext)- At the Start of the pulse 
(t=0). Yd increases, following a canonical Resistance- 
Capacitance (RC) charging curve. Once Yd surpasses a 
threshold voltage (which we will call Vd-.imt)^ it tran- 
sitions sharply from increasing Vd to decreasing. We 
assign this change to an IMT event occurring in the 
VO2, which effectively lowers the internal shunt resis- 
tance K]j of the capacitor. With lower internal resis- 
tance, the capacitor undergoes rapid discharge and Vd 
plummets. This discharge continues until V^ reaches a 
lower threshold voltage {Vd-.mit at which a second event 
- which we similarly assign to a a Metal-to-Insulator- 
Transition (MIT) - restores the high internal device re- 
sistance. The process reverses and this sequence of events 
repeats; alternating charging and discharging between 
IMT and MIT events with a fairly stable periodicit}^^. 



III. GRAIN-SCALE MODEL OF OSCILLATIONS. 

Our hope is that by developing a model for these ob- 
served oscillations, we may gain insight into the driving 
mechanism behind them. We start with a 2D network 
of square VO2 grains of differing size, each with total 
resistivity dependent on its size and its state (metal or 
insulator) T he g ranularity of polycrystaline VO2 

is well documente d^ ^ although the size of grains may 
vary considerably from one VO2 preparation to another. 
There is also evidence to suggest the percolation length 
scales for the IMT may not always coincide with the crys- 
tal granularit}'^^ A cartoon illustrating our model ar- 
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FIG. 3. Cartoon of the grains resistor-network. Left and 
right sides are the beginning of the nickel electrodes. In the 
cartoon, the size of squares represent a distribution of grain 
sizes and resistor colors serve to highlight the corresponding 
different resistance. 



rangement is shown in Figure |3] The network consists of 
N^xNj grains, and in our model we restrict our investiga- 
tion to a 50x50 network array to keep computation time 
manageable. 

This grain network is placed in an external circuit con- 
taining resistance Rext and capacitance Cg^t, and driven 
by Vapp, as shown in Figure [T] The circuit differential 
equation (Equation [T]), is solved via Runge-Kutta time- 
stepping. 
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At each time- step, we solve for the internal state of the 
grain network. This includes solving a Kirchoff net- 
work problem^^ for the voltage across each grain and the 
Thevenin effective circuit resistance Rd- The other VO2 
effective circuit parameter Ci:) is found by a differential 
capacitance equation (Equation [2| which can be evalu- 
ated via a self-consiste nt B ruggeman effective medium 
formulation (Equations |3|4|) 
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Superscript i and j are row and column indices for the 
grains (running to and Nj total). The binary ma- 
trix X*-^ = 1 if the grain i,j is metal and if insulator. 
?7 is a capacitive fractional-fields factor (the proportion 



of the device capacitance which involves the VO2 dielec- 
tric), which is found via finite element simulation using 
the COMSOL commercial code package. Co is a geomet- 
rical capacitance which is determined empirically, fitting 
1/RC to the capacitive charging curve. External circuit 
parameters such as Rext and Vapp are taken directly from 
experimental values. The extrema values for R/) (metal) 
and Ri:) (insulator) are taken from temperature data. 



A. Thermal triggering 

The temperature driven IMT-MIT has been investi- 
gated in great detail and we ground our model using 
experimental data giving resistance as a function of tem- 
perature R(T) through the phase transition. This data is 
shown in Figure |4^, and displays the characteristic sharp 
change in resistivity around 345K. 

T6l26l 



Then, in a procedure similar to previous works ' , we 
assume each grain will undergo an IMT in response a 
"high-threshold" temperature T)^^^, and a MIT at a 
low-threshold T^j^. 
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We assign a stochastic distribution to the values of Tj^j^j. 
and T^jj. throughout the network, following the Gaus- 
sian form 
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The values of Tq and are fit to the experimental R(T) 
data shown in Figure H^. From this, we find a variance 
of To = 340^K and cr|^. l*To reproduces the shape of 
the IMT fairly well. This fit gives us confidence that we 
understand the thermal response of VO2, even though 
our thermal model is simple compared to some previous 
treatments. 

Using this temperature-only-triggered model, we at- 
tempt to reproduce the oscillations shown in Figure [2j 
Our model tracks the power dissipated in each grain and 
employs a finite-element method to solve for the grain 
and substrate temperatures as a function of time. Ma- 
terial thermal parameters are taken from literature, and 
the enthalpy of phase-transition for VO2 is included. Fig- 
ure |4]3 re-plots our experimental Vd oscillations in blue, 
and the model results in green. The result is striking: al- 
though we can track Vd for a short while, we are unable 
to produce any oscillatory phenomena. 

Generalizing the behavior that prohibits oscillations, 
thermal-initiated IMTs exhibit a runaway behavior 
rather than the self-stabilizing oscillatory nature seen in 
Figure [2] This is apparent in the average VO2 grain 
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FIG. 4. a) Experimental resistance of our device as a function 
of temperature Rd(T) (black). Ri:)(T) results from our model 
are overlaid for a best-fit value of variance: cf^mt — 0-1 * 
(see Equation |7| , which replicates the observed temperature- 
driven IMT fairly well. Also shown is a poor-fit result for 
c^iMT — 0-01 * ^0, given to show the effect of this parameter, 
b) Attempt to replicate oscillations using thermal-only trig- 
gering, with values of PiTy^rp) from above. We plot Vd from 
our model (green) overlaid on experimental Vd (blue). We 
observe no oscillations, only a single IMT event. We also plot 
the average VO2 temperature from our model, illustrating 
a runaway heating behavior that precludes oscillations. 



temperature plot (red) of|4]3. Note that although indi- 
vidual different grains may attain different temperatures 
over the course of oscillations, such gradients equalize 
quickly within the network. Average temperature re- 
mains a fairly accurate and easily visualized metric of 
the VO2 oscillation thermodynamics. As the thermally- 
triggered IMT occurs, VO2 temperature skyrockets even 
while Yd discharges. Once the temperature of the dis- 
charged device has surpassed Tq^^^, there is no way for 
it to cool thereafter. Power dissipation and temperature 
both settle towards steady-state with the device firmly 
in the metallic state. The causes of this process will be- 
come more apparent as we discuss thermal dynamics in 
Section [mDl 
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FIG. 5. The experimental data from Figure |2]replotted (blue) 
with numerical results from our model overlaid (green). This 
electrical-and-thermal triggering model replicates the experi- 
mental data quite well, tracking the oscillation periodicity and 
producing similar transitions at both Vd-.imt and Yd-.mit- ) 



tribution as we did for ^^jmt ^^mit- 
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We lack direct data to which to fit these distributions 
(as we did in Figure |4|i). Thus, we retain the same value 
for the variance found in above, <jy=0.1*Vo. From an 
energetics perspective, this makes a great deal of sense 
- both distributions are surely tied to the same underly- 
ing Mott-physics. Vq remains a fitting parameter in our 
model. 

As we have not removed possible thermal triggering, 
the conditions for grain transition can now be stated as: 
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This combined triggering criterion reproduces both the 
waveshape and periodicity of the experimentally ob- 
served oscillations quite well. In Figure [5| we re-plot the 
experimental data from Figure |2] in blue, and numerical 
results from our model are shown overlaid in green. 



B. Electrical triggering 



C. Voltage-temperature dependence 



With the failure of a temperature- on/^-triggering 
model to produce oscillations, and following insights 
from previous work on voltage-induced effects in 
YQ^ |22i2 6 | 37 | 38| introduce an electric-field driven 

transition. To do this, we also assign a voltage drop at 
which each grain undergoes phase transition Y^j^t 
^^MiT' again use random values from a normal dis- 



Although section [ni B| demonstrated voltage is the pri- 
mary trigger, device temperature still plays a role in os- 
cillations. There is a known dependence of oscillation 
amplitude on device temperature^^. If we look carefully 
at the data in Figure [5) we notice a subtle decay enve- 
lope to the amplitude of the Vjmt oscillation peaks. We 
believe this envelope is caused by a thermalization of the 
device on a multi-oscillation timescale. To accommodate 
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FIG. 6. Experimental data (blue) for oscillations over a long 
100/is Yapp pulse. We observe a clear decay envelope to the 
amplitude of V/mt in the oscillations. For clarity, the bottom 
half of the oscillations are omitted as we observe no enveloping 
here. To show this decay is a thermalization envelope caused 
by device heating, we also plot the temperature dependent 
Voj^^(T) from equation 13 (red). 



this, we include a temperature-dependence to the IMT 
transition voltage as 



ViiiAT) ^ {'^{T - To) 



IMT 
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where To=295^K. is a liner temperature coefficient, the 
fitting of which we discuss below. Without including this 
voltage-temperature interplay, our model quickly loses 
sync with the experimental data over the course of several 
oscillation periods. This thermalization envelope is most 
clearly observed over a long pluse, and experimental data 
(solid blue) for a 100/is pulse is shown in Figure |6] For 
clarity in this figure, the bottom half of the oscillations is 
omitted from view (interestingly, we observe no envelope 
of Vm/t). 

Using the data from Figure |6] combined with our 
thermal- finite element model, we can fit a value for k,. 
Overlaid on the experimental data (dashed line) is this 
fit, which gives = —8.3x10"^. This value of is used 
in our oscillation model. 

Explanation of this a long time-scale thermalization is 
straightforward. Although locally the VO2 film may heat 
or cool quite quickly in response to current through its 
volume, the 500/im thick sapphire substrate is compara- 
tively massive. The large thermal inertia of the substrate 
smooths out oscillatory heating away from the film; and 
when heated only from the top the substrate can require 
tens of micro-seconds to reach a steady-state temperature 
gradient. 



D. Temperature dynamics 
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FIG. 7. Model results of temperature dynamics during \ d os- 
cillations. The horizontal dashed line at T/mt=340^K shows 
mean transition temperature, and demonstrates that the av- 
erage temperature reached during oscillations falls short of 
what's needed to drive the IMT transition. We also divide 
the oscillation into insulating/charging regions(I) and metal- 
lic/discharge regions (II). This draws attention to the sharp 
increase in temperature which occurs only ajier the Yimt 
trigger, qualitatively discrediting a thermal-driving picture 
for oscillations. 



To begin, in Figure [7| we plot the average grain tem- 
perature (red) during oscillations along with Vd (green) 
- both from our model. Looking at Figure [7| quantita- 
tively, we notice that the average temperature reached is 
not sufficient to trigger oscillations. Although the peak 
average temperature during the discharge cycle of the os- 
cillation comes close to reaching T/mt, a thermal-driving 
event would have to occur at-or-before the peak in V^. 
The model results indicate that the device is well be- 
low T/MT when Yd-.imt is reached. The temperature 
range in Figure [7| agrees fairly well with previous nu- 
merical work on the subject^, and our own investiga- 
tions using commercial finite element package COMSOL. 
Compared to COMSOL, our home-grown finite-element 
code over-estimates temperatures reached - perhaps due 
to difficulty in modeling all of the spatially massive sub- 
strate (this is not an issue for commercial packages such 
as COMSOL). However, commercial codes cannot be run 
from within Runge-Kutta time-stepping of equation [l] 
and our code thus critically allows us to calculate tem- 
perature dynamics during oscillations. 

Although a quantitative argument for non-thermal 
triggering seems compelling, we are acutely aware that 
precisely solving for temperature can be difficult in such 
nanoscale systems. Material properties can differ from 
published bulk values, and interface effects can dominate 
transport and heatin^^. As we look closely though, the 
power dissipation in the device also appears qualitatively 
unfit to explain the oscillations. A thermally-driven tran- 
sition would have to follow the logical sequence: 



In the above Section llll Cl we have discussed the critical 
interplay between temperature and voltage-triggering, 
showing how a long multi-oscillation timescale thermal- 
ization envelopes the oscillation amplitude. In this sec- 
tion, we look closer at the temperature evolution on the 
timescale of the oscillation period. 



I The insulating device heats with applied Vjj until it 
reaches T/mt, where it undergoes an IMT (becom- 
ing metallic). 

II The metallic device discharges its stored capacitive 
energy through the its own volume - cooling as it 
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discharges - until it reaches Tmit where it undergoes 
MIT (becoming insulating). 

Ill The process repeats. 

However, Figure [7| illustrates that region II which oc- 
curs after Yd-.imt is a region of maximum power dissipa- 
tion; a region of heating not cooling. A simple ohms-law 
argument explains: Just before and just after the IMT, 
Vi) is approximately Yd-.imt- However, the resistance 
has changed by a factor of 10, and thus the power 
dissipated (P=V|).j^^/Ri:)) is substantially greater dur- 
ing region II than during region I. Thus, we come to 
the conclusion that a purely-thermal explanation for the 
oscillations is qualitatively as well as quantitatively mis- 
matched to experimental data. 



Summarizing Section HI, we have identified voltage 
as a key player in triggering observed oscillations on 
the grounds of several thermal arguments. One in- 
teresting question then is whether electrostatic voltage 
may also trigger the also IMT in a current-free (ie. 
FET) configuration. The joule-heating present in our 
two-terminal device complicates matters, in light of the 



voltage/temperature interplay identified in Section IIIC 



Previous work h as suggested such electrostatic switching 
can exist although these early results await further 

confirmation. Exploring the phase-space defined by the 
interplay of temperature, electrostatic field, and current 
in these VO2 oscillations may reveal information about 
the correlated electron dynamics and energy scales asso- 
ciated with the Mott-transition. 



IV. PERCOLATION 

In this section, we go into further detail on the mech- 
anisms of the IMT and MIT transitions. Polycrystalline 
VO2 is known to exhib it percolative behavior during 
phase transition^^^EHlll^ and this has interesting effects 
on a voltage-triggered transition. Using the model from 



Section [Hi] which accurately predicts the observed elec- 
trical oscillations, we attempt to gain insight on several 
of the internal processes during oscillatory events. 



stochastic distribution for Y^j^t (^^^ Equation M), this 
can be a random grain anywhere in the network. In ex- 
periments, it is seen that the phase-transition is often 
seeded at particular places such as defects or boundaries. 

The unlucky grain that first hits its IMT trigger condi- 
tion undergoes an IMT. Once this grain becomes metal- 
lic, it supports a lower voltage drop {Kins/^met ^ 20) 
- which shifts much of its voltage burden to neighbor- 
ing grains. The neighboring grains in turn become in- 
creasingly likely to undergo their own IMT events. The 
IMT spreads across the entire sample in an avalanche- 
like manner. This process is depicted in Figure for a 
network of 50x50 grains. The upper sequence of black 
and white frames shows whether each grain is insulat- 
ing (white) or metallic (black). The lower color frames 
depict the voltage drop across each grain. The neighbor- 
neighbor grain interaction, mediated by voltage drop is 
quite evident. 

As is common in percolative systems, the Thevenin 
resistance oi the network is quite sensitive to the 
spatial distribution of triggered grains. The device resis- 
tivity Ri:) is shown above each frame in the sequence of 
Figure [8) and we see the largest drop occurs in frames #6 
and 7, where the percolation path is completed from left 
to right. Once a conducting path forms Rd plummets, 
and Vd begins to drop as the device discharges. 

During discharge, there comes a point where drops 
far enough that a similar process happens in reverse. This 
MIT is depicted in Figure (8)3. However, the MIT process 
is not exactly the reverse of the IMT. The equations gov- 
erning resistors in parallel tend to "favor" low-resistances 
in the following manner: Decreasing the value of one re- 
sistor (in a parallel network) lowers the Thevenin resis- 
tance significantly, but raising the value of a single re- 
sistor has only a little effect on the Thevenin resistance. 
For this reason, the avalanche-like behavior observed in 
Figure [8^ is not seen in Figure [8]3. Instead, the process 
much more closely resembles random percolation, with 
only moderate neighbor-neighbor interaction. This dif- 
ference in the mechanisms between IMT and MIT may 
explain the difference in sharpness of the transitions at 
Vdimt ^DiMT observed in experiments.'^ 



A. Percolative avalanche driven oscillations 

In several previous work4^^^, avalanche-like MIT and 
IMT transitions have been observed under the right con- 
ditions. The immediacy of the observed change from 
charging to discharging in our oscillations leads us to sus- 
pect similar avalanche behavior in our electrically-driven 
device. Examining the details of our model, we see that 
the voltage drop across any grain in the network (see Fig- 
ure [3| is proportional to the resistance of the grain. Dur- 
ing the charging cycle of the waveform voltage across the 
entire device (V^) increases, and V*-^ across each grain 
does as well. This charging continues until one "unlucky" 
grain hits its Y^j^t fii"st. Because the grains receive a 



B. Effective medium effects 

The percolative nature of the VO2 phase transition al- 
lows for an inhomogeneous intermediate state where both 
metallic and insulating VO2 coexist. The dielectric con- 
stants of metal and insulating phase VO2 are distinct; 
and when both phases can be present in a composite it 
leads to interesting properties. The average response of 
the inhomogeneous sample is described by an effective 
medium, and can have radically different values than ei- 
ther. This leads to quite interesting and novel effects. For 
example, the inhomogeneit}^^^ of polycrystalline VO2 
mid-transition is responsible for observed memristanc^^^ 
and memory-capacitance^^ 
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FIG. 8. Step-by-step depiction of the avalanche- like transition for a 50x50 grain network. The time of each frame increases 
from left to right, a) shows the IMT transition occurring at Vd-.imt, giving a bi-color plot (top) indicating whether each grain 
is metallic or insulating and the the voltage (bottom) across each grain in the network, b) shows the same plots for the MIT 
transition occurring at Vd-.mit- 



The same memory-capacitance as reported iiP^l has 
previously been attributed as playing a key role in 
voltage-controlled oscillations in VOj^. This is a ques- 
tion we are situated to investigate in more depth using 
our model. To look closely at the effects of capacitance 
on the oscillations, we focus attention on a single IMT 
transition event. In Figure [9j we plot Cd{^) (as calcu- 
lated from equations [2][4|, along with the familiar 
and Rd- 

Looking at Figure [9) as VO2 transitions from insulat- 
ing to metallic at the IMT, Kjj drops monotonically to 
its metallic-state value. The capacitance Cd, however, 
briefly increases before also decreasing to its metallic- 
state value. This increase is due to the coexistence of 
metallic and insulating grains, and is predicted by effec- 
tive medium (Equation [4]) . However, as Figure [9] shows, 
the increase in is a small effect, and is contained to a 
short timespan near the start of the IMT. This leads us 
to believe the effective medium behavior of Cd has only 
a minor influence on the shape of the oscillations, and is 
not a primary driver. We observe the same increase in 
Cd Sit the MIT transition edge, but it also is too small 
and short-lived an effect to bear responsibility for the 
oscillations. 
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FIG. 9. Model values for Cd plotted with Vd and Kd over 
one oscillation period. Effective medium within the VO2 (see 
Equation [2]) causes a brief spike in. Cd at immediately near 
the IMT and MIT transitions. As seen, this spike effect is 
quite small compared to the overall change of both Cd and 
Rd, and exists only for a very brief fraction of the oscillation 
period. The jagged shape of the model curves for Rd and 
Cd reveal the small jumps typical of percolation discussed in 
Section [IVA] 



V. SUMMARY 



In this work, we have discussed the thermal and electri- 
cal driving mechanisms behind observed oscillations oc- 
curring in VO2 films. In addition to experimentally con- 
firming the oscillations reported by Y.W. Lee et.al.^, we 
have compiled a numerical model which is able to repli- 
cate and explain these oscillations in terms of an voltage- 
triggered Insulator-to-Metal phase-transition. Tempera- 
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ture is known to trigger the IMT in VO2, and tempera- 
ture plays some role in the shape of the oscillations. How- 
ever, we find that temperature-only triggering cannot ex- 
plain oscillations. This result is likely of great importance 
for applications of these oscillations, as repeated thermal- 
cycling typically appreciably shortens device lifetime. 

One question that remains unaddressed is the role, if 
any, of the structural phase transition in an electric-field 
triggered transition. The temperature driven IMT in 
VO2 exhibits a structural transition that happens concur- 
rent with the electronic reconfiguration. However, there 
is evidenc d^'^'^^ ' ^^^^^' to suggest the electronic and struc- 
tural transitions are not necessarily linked, but merely 
overlaid. This structural electronic decoupling suggests 
the electronic correlations in VO2 play an important role 
in the IMT phase transition. 

Studies which probe the crystal structure simultane- 
ous with these oscillations have not yet been reported 
- likely because the time and length scales associated 
with the VO2 oscillator devices greatly complicate exper- 
iments such as x-ray diffraction. However, as mentioned 
in the introduction (Section |l| the question of whether 
or not the structural transition occurs has great impli- 
cations about the longevity of these devices. In many 
applications, devices could easily be expected to perform 
10^^ to 10^^ oscillation events over their lifetime, a likely 
impossibility if crystallographic changes are occurring. 
An evident goal for the new future is to experimentally 
investigate the existence of structural transition in oscil- 
lations in the near future. 
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